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1. Introduction 


The fractional calculus of variations deals with optimization of functionals that 
depend on some fractional operator Hi- This is a fast growing subject, and 
different approaches have been developed by considering different types of La- 
grangians, e.g., depending on Riemann-Liouville or Caputo fractional derivatives, 
fractional integrals, and mixed integer-fractional order operators @-0j. The com¬ 
mon procedure to address such fractional variational problems consists in solving 
a fractional differential equation, called the Euler-Lagrange equation, which every 
minimizer/maximizer of the functional must satisfy. With the help of the boundary 
conditions imposed on the problem at the initial time t = a and at the terminal 
time t = b, one solves, often with the help of some numerical procedure, the frac¬ 
tional differential equation and obtain the possible solution to the problem |10H15|]. 
When one of the boundary conditions is missing, that is, when the set of admissible 
functions may take any value at one of the boundaries, then an auxiliary condition, 
known as a transversality condition, needs to be obtained in order to solve the frac¬ 
tional equation [T^]. In this paper we do not only assume that x(b) is free, but the 
endpoint b is also variable. Thus, we are not only interested in finding an optimal 
curve x(-) but also the endpoint of the variational integral, denoted in the sequel by 
T. This is known in the literature as a free-time problem, and is already well stud¬ 
ied for integer-order derivatives (see, e.g., 0 )- Moreover, we consider a combined 
Caputo fractional derivative, that involves the left and the right Caputo fractional 
derivative [H-0. Also, the order of the fractional operator is not a constant, and 
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may depend on time 21 [231]. Very useful physical applications have given birth to 
the variable order fractional calculus, for example in modeling mechanical behav- 
23-0. adays, variable order fractional calculus is particularly recognized 


lors 


as a useful and promising approach in the modelling of diffusion processes, in order 
to characterize time-dependent or concentration-dependent anomalous diffusion, or 
diffusion processes in inhomogeneous porous media [13]. The importance of such 
variable order operators is due to the fact that the behaviors of dynamical systems 
often change with their own evolution [261 ]. For numerical modeling of time frac¬ 
tional diffusion equations we refer the reader to 245 • Our variational problem may 
describe better some nonconservative physical phenomena, since fractional opera¬ 
tors are nonlocal and contain memory of the process 


27-29]. We believe that it is 


reasonable that the order of the fractional operators is influenced also by the acting 
dynamics, which varies with time, and we trust that the results now obtained can 
be more effective to provide a mathematical framework to some complex dynamical 
problems [US]. 

The paper is organized in the following way. In Section [2] we introduce the basic 
concepts of the combined variable-order fractional calculus, and we recall two for¬ 
mulas of fractional integration by parts (Theorem [2T]). The problem is then stated 
in Section [3l consisting of the variational functional 


J{x, T) = f L (t, x(t), c dt + 4>(T, x(T)), 


where c ’ ^x(t) stands for the combined Caputo fractional derivative of 

variable fractional order (Definition [lj) , subject to the boundary condition x(a) = 
x a . Our goal is to find necessary optimality conditions that every extremizer (x, T ) 
must satisfy. The main results of the paper provide necessary optimality conditions 
of Euler-Lagrange type, described by fractional differential equations of variable 
order, and different transversality optimality conditions (see Theorems 13.11 and 
13.21) . Some particular cases of interest are considered in Section [H We end with 
two illustrative examples (Section [5]). 


2. Fractional calculus of variable order 


In this section we present the fundamental notions of the fractional calculus of 
variable order. We consider the order of the derivative and of the integral to be a 
continuous function «(•,•) with domain [a,6] 2 , taking values on the open interval 
(0,1). Following [U], given a function x : [a,b] -^Mwe define: 

• the left and right Riemann-Liouville fractional integrals of order «(•,•) by 


and 


J^xit) = f 1 (t - rr^-'x(T)dT 
J a r(a(i, t)) 


t-£ (v) ®(t)= [ p7 T TyT ( r ~ l x{r)dT , 

Jt r(a(TD) 


respectively; 
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the left and right Riemann-Liouville fractional derivatives of order a(-, •) by 


D? W) x{t) = ~T7 


d 


dt J a r(l - a(t, T )) 


(t - T)~ a ^x{r)dT 


and 




-1 


t T(1 -a{r,t)) 


(r - t)~ a ^x(T)dT, 


respectively; 

the left and right Caputo fractional derivatives of order a(-, •) by 




1 


I a r(l - a(t, T )) 


(t — t) a ^’ T ' ) X^ (r)dT 


and 


?Df^x(t) = 


-1 


<t r(l -a(r,t)) 


(r-i)- Q ( T ’ t )i( 1 )(r)rfr, 


respectively. 

Motivated by the combined fractional Caputo and Riemann-Liouville definitions 
IH0, we propose the following definitions. 

Definition 1 Let a, ft : [a,b ] 2 —>• (0,1) and 7 = ( 71 , 72 ) £ [0, l] 2 . The combined 
Riemann-Liouville fractional derivative operator Illy ,) - g ( j e fj ne( j by 

= 7 lajD «(v) +72tjD ^v) ; 
acting on x £ C([a, b\) in the following way: 

= 7l a H“ ( '’' ) x(t) +72 tD^'^xit). 

Analogously, the combined Caputo fractional derivative operator, denoted by 
c is defined by 

= 7l C D a(;-) + 72 C^(v) ; 

acting on x £ C 1 ([a, b]) as expected: 

C D<DM-D x (t) = ^D^ W) x{t) +72?D^’- ) x{t). 


The following theorem is proved in [35] and is a generalization of the standard 
fractional formula of integration by parts for a constant a (see, e.g., formula (2) of 

0 ). 
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Theorem 2.1 (Theorem 3.2 of (35|) //i,y£ C^fa, b\, then the fractional integra¬ 
tion by parts formulas 


y( t )a D t < " )x ( t ) dt = [ x{t) t Dl { '''\{t)dt+\x(t) t ll " (v) y(i) 


t=b 


t=a 


and 



?x(t)dt 


x(t) a D^'''\{t)dt 


x(t) a ll a(v) y(t) 


t=b 


t=a 


hold. 


3. Necessary optimality conditions 


Consider the norm defined on the linear space C 1 ([a, b \) x M by 


ll(M)ll 


max \x(t) \ + max 

a<t<b a<t<b 


c D^^x(t) 


+ \t 


Let D denote the subset C ,1 ([a, b ]) x [a, b\ endowed with the norm ||(-, -)||, such that 
c exists and is continuous on the interval [a, b\. 

Definition 2 We say that (x*,T*) £ D is a local minimizer to the functional J : 
D —>• M if there exists some e > 0 such that 


V(s,T) G D : ||(s*,r*) - (x,T )|| < e =7 J(x*,T*) < J(x,T). 

Along the work, we denote by diZ, i £ {1,2,3}, the partial derivative of a func¬ 
tion z : M 3 —> M with respect to its zt.li argument, and by L the Lagrangian 
L : C 1 ([a, b] x M 2 ) —> M. For simplicity of notation, we introduce the operator 

defined by 


[x}^(t) = (t,x(t), c D^^x(t)) ■ 

Consider the following problem of the calculus of variations: hnd the local mini- 
mizers of the functional J : D —> R, with 

J{x, T) = [ T L[x]^(t)dt + x(T)), (1) 

J a 

over all (x, T) £ D satisfying the boundary condition x(a) = x a , for a fixed x a £ M. 
The terminal time T and terminal state x(T) are free. The terminal cost function 
4> : [a, b] x R — > M is at least of class C 1 . 

The next theorem gives fractional necessary optimality conditions to problem 
(fTj) . In the sequel, we need the auxiliary notation of the dual fractional derivative: 

D P(;-),<*(;■) = ry 2a £)hi;-) + 7l where 7 = ( 72 , 71 )- 
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Theorem 3.1 Suppose that (x,T) is a local minimizer to the functional (JT]) on D. 
Then, (x, T ) satisfies the fractional Euler-Lagrange equations 


d 2 L[x]^(t) + D^’^dsLix^it) = 0, (2) 

on the interval [a,T], and 

72 (aD^dzL[x\^{t) - T D^d 3 L[x] a /{t)) = 0 , ( 3 ) 

on the interval [T,b ]. Moreover, (. x,T ) satisfies the transversality conditions 


L[x\^(T) + d l 0(T, x(T)) + d 2 fi{T, x(T))x'(T) = 0, 
7 i d 3 L[x]^(t) - 72 v d 3 L[x]°’P(t) 

72 d 3 L[x\^ {t) - d 3 L[x]j ,f3 (t) 


J t=T 
= 0 . 


+ d 2 ^(T,x(T)) = 0, 


- t=b 


(4) 


Proof. Let (x, T ) be a solution to the problem and (x + eh, T + eAT) be an admis¬ 
sible variation, where h £ C fl ([a, 6 ]) is a perturbing curve, AT £ M represents an 
arbitrarily chosen small change in T and e £ M represents a small number (e —> 0). 
The constraint x(a) = x a implies that all admissible variations must fulfill the 
condition h{a) = 0. Define j(-) on a neighborhood of zero by 


j (e) = J(x + eh,T + eAT) 

/*T+eAT 

= / L[x + eh]% fi (t)dt + <j>(T + eAT,(x + eh)(T + eAT)). 


The derivative j'(e) is 

/*T+eAT . v 

/(e) = J (3 2 T[x + eh]“ ,/3 (t)h(t) + <9 3 T[x + eh}^{t) c dt 

+ L[x + eh]°’P(T + eAT)AT + (T + eAT, (x + eh){T + eAT)) AT 
+ d 2 fi (T + eAT, (x + eh)(T + eAT)) (x + e/i)'(T + eAT). 


Considering the differentiability properties of j, a necessary condition for (x,T) to 
be a local extremizer is given by j'(e) | e=0 = 0 , that is, 

j^ (d 2 L[x]“’P(t)h(t) +d 3 L[x]^(tfD^’^h(t)) dt + L[x]«’P(T)AT 

+ d 1( f (T, x(T)) AT + d 2 (/)(T, x(T)) [/r(t) + x'(T) AT] = 0. (5) 

The second addend of the integral function ©, 

d 3 L[x]^(t) c D^’^h(t)dt, ( 6 ) 
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can be written, using the definition of combined Caputo fractional derivative, as 


d 3 L[x]^(t) c D^’^h(t)dt 


f d 3 L[x]^(t) [ 71 + 72 

J a L 


dt 


= 71 / d 3 L[x]^(Z)fZ}“ (v) /i(Z)df 


+ 72 


[ d 3 L[x]^(t)f Dl w) h(t)dt - f 9 3 L[x]"’ /3 (f) f c D% ( ' r) h(t)dt 

J a J T 


Integrating by parts (see Theorem 12.11) . and since h(a) = 0, the term <(U]) can be 
written as 


7i 


J t=T 


+ 72 


h(t) t D^d 3 L[x]*’P(t)dt + [h(t) t 4- a ^d 3 L[x]*’P(t) 

[ h(t) a D^' , ' ) d 3 L[x]^(t)dt - /i(t) Q / t 1_/3( '’' ) 5 3 L[x]"’ /3 (t) 

J a 

- (^l\(t) T D^d 3 L[x}^(t)dt - \h(t) T ll-^’- ] d 3 L[x]«’P(t) 

h(t) T ll-^d 3 L[x]^(t) 


t=b 


J t=T 


Unfolding these integrals, and considering the fractional operator Z+v ’ l ' rr ' ’ ^ with 
7 = (72,7i)) then (|6|) is equivalent to 


[ T h^D^^dsLixj^^dt 
J a 


f 72 h(t) [ a D^d 3 L[x]^(t) - T D^d 3 L[x]^(t) 
71 t I l T -^^d 3 L[x\^(t) - 72 T ll~ PW) dzL[x]^ 
T ll-^d 3 L[x}^(t) - 


dt 


J t=T 


- t=b 
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Substituting these relations into equation (f5|). we obtain 


0= [ T h(t ) \d 2 L[x]^(t) + 
Ja *■ 


dt 


+ 


I2h{t) \ a D^d 3 L[x\^{t) -TD^dsLix^it) 


dt 


+ h(T ) 71 Jt a(v) d 3 £M“ ,/3 (£) ~l 2 Tl] d 3 L[x]^ (t) + d 2 (j)(t,x{t)) 


J i=T 


+ AT 


T[x]" ,/ 3 (t) + x(t)) + d 2 cj)(t, x{t))x'(t) 


+ h(b) [72 (Tl^^dsLixp^t) - alt-K^dsLlxV’P 


J t=T 

a 
7 


t=b 


(7) 


As h and AT are arbitrary, we can choose AT = 0 and h(t ) = 0, for all t G [T, 6], 
but h is arbitrary in f G [a, T). Then, for all t G [a, T], we obtain the first necessary 
condition (|2|) : 


d 2 L[x]«’P(t) + D^’- ) ’ a ^ ) d 3 L[x}^(t) = 0. 

Analogously, considering AT = 0, h(t) = 0, for all t G [a, T], and h arbitrary on 
(T, b], we obtain the second necessary condition ©: 

72 (aDf M d3L[x]^(t) -TD% W) d3L[x]°’P(t)) =0. 

As (x,T) is a solution to the necessary conditions © and ©, then equation (J7]) 
takes the form 


0 =h(T) 7 ! “ ( ' , ' ) 5 3 T[x]“’ / 3 (f) - 72T7 * 1 / 3 ( ' , ' ) «9 3 -b[x]“’ / 3 (f)+5 2 ^(t,x(t)) 


J t=T 


+ AT 


T[x]" ,/3 (t) + d\(f>(t, x{t)) + d 2 (f>(t , x(t))x(t) 


J t=r 


+ M&) [72 (tI] /3(v) d 3 L[z]“’ /3 (f) - all ^dzLlx 


- t=b 


( 8 ) 


Transversality conditions © are obtained for appropriate choices of variations. □ 


In the next theorem, considering the same problem (jTJ) , we rewrite the transver¬ 
sality conditions © in terms of the increment on time AT and on the consequent 
increment on x and A xt, given by 

Ax t = (x + h)(T + AT) -x(T). (9) 

Theorem 3.2 Let (x,T) be a local minimizer to the functional © on D. Then, 
the fractional Euler-Lagrange equations (ED and © are satisfied together with the 
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following transversality conditions: 


(L[xp P (T) + dKl>(T,x(T)) 

+x'(T) ’dsL[x]“’P(t) - v d 3 L[x]j’^ (t) 


7i tlx a<K '^d 3 L[x]^(t) - 72 T l] /3C "’" ) «9 3 L[x]^(t) 
T / t 1 _/ 3 (v) d 3 L[xp f} (t) - a / t 1 _/ 3 (v) 9 3 L[x]“ ,/ 3 (t) 


= 0 , 


t=T 


72 


t=T 

= 0 . 


+ <9 2 <KT,x(T)) = 0 , 


- t=b 


( 10 ) 


Proof. The Euler-Lagrange equations are deduced following similar arguments as 
the ones presented in Theorem 13.11 We now focus our attention on the proof of 
the transversality conditions. Using Taylor’s expansion up to first order for a small 
AT, and restricting the set of variations to those for which h'(T ) = 0, we obtain 


(x + h)(T + AT) = (x + h)(T) + x'(T)AT + O(AT) 2 . 


Rearranging the relation ([9]) allows us to express h(T) in terms of AT and Axt- 

h(T) = Ax T ~ x'(T)AT + O(AT) 2 . 


Substitution of this expression into (JH]) gives us 


0 =A x T 7 i i /^“ ( '’' ) 0 3 L[x]"’ / 3 (t) - l2T ll~ P{w) d^L[x\^{t) + d 2 (j){t, x(t)) 
+ AT T[x]“ ,/ 3 (t) + di4>(t, x(t)) 

-x'(t) (71 t /^“ (v) d 3 T[x]“’ / 3 (t) - v d^Llx]^ 


J t=T 


J t=T 


+ h(b ) [72 (T/ i 1_/3(v) 9 3 L[x]^(t) - a ll-^- ) d 3 L[x}^(t))] + O(AT) 2 . 


- t=b 


Transversality conditions (HOD are obtained using appropriate choices of variations. 

□ 


4. Particular cases 

Now, we specify our results to three particular cases of variable terminal points. 


4.1. Vertical terminal line 

This case involves a fixed upper bound T. Thus, AT = 0 and, consequently, the 
second term in (JH]) drops out. Since Axt is arbitrary, we obtain the following 
transversality conditions: if T < b, then 


7i “ (v) d 3 £M7 ,/3 (£) - l2Tl] /3( ' , ' ) <9 3 T[x]“’ ,3 (t) 
<9 3 L[z]“ ,/ 3 (t) - a / f 1 ' / 3 ( ' , ' ) 9 3 L[x]"’ / 3 (t) 


t=b 


t=T 

= 0 


+ d2<j>(T, X (T)) = 0, 
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if T = b, then A xt = h(b) and the transversality conditions reduce to 


7i tJr " (v) d 3 T[x]"’ / 3 (t) - 72 tl? / 3 ( ' , ' ) <9 3 T[x]"’ / 3 (f) + d 2 (t>{b, x(b)) = 0 


J t=b 


4.2. Horizontal terminal line 

In this situation, we have A xt = 0 but AT is arbitrary. Thus, the transversality 
conditions are 


'L[x\^{T) + d l( l>{T, X (T)) 

+x'(T) [ 72 T 1 7 f 1 ~ / 3 ( '’' ) <9 3 T[x]"’ / 3 (t) - 7 it/^ Q( '’' ) 9 3 T[x]“’ / 3 (t) 


72 


t/j 1 / 3 ( ' , ' ) <9 3 T[x]"’ / 3 (f) - j] / 3 ( ' , ' ) <9 3 L[x]"’ / 3 (f) 


= 0 , 


J t=T 


= 0 . 


J t=b 


4.3. Terminal curve 

Now the terminal point is described by a given curve 1 p : C fl ([a, 6 ]) —> M, in the 
sense that x(T) = From Taylor’s formula, for a small arbitrary AT, one has 

A x(T) = ip'(T)AT + O(AT) 2 . 


Hence, the transversality conditions are presented in the form 


L[x}^(T) + 3k j>(T, x(T)) + d 2 <P(T, x(T)W(T) 

+ {x'(T) - ip'(T)) [72 ’ } 9 3 T[x]“’ / 5 (f) - 71 i Iy'" ( ' , ' ) 3 3 L[x]"’ / 3 (f) 


72 


T l]~ Kw) d^L[x]^{t) - a 7 t 1 - /3(v) 3 3 T[x]“’ /3 (f) 


J t=T 


= 0 , 


= 0 . 


J t=b 


5. Examples 

In this section we show two examples for the main result of the paper. Let a(t, r) = 
a(t) and /3(t,r) = f3(r ) be two functions depending on a variable t and r only, 
respectively. Consider the following fractional variational problem: to minimize the 
functional 


J(x,T) 



2 a(t) — 1 + 



jX—a(t) 

2r(2 -a(t)) 


(10 

2T{2-m) 


21 


dt 


for t £ [ 0 , 10 ], subject to the initial condition x( 0 ) = 0 and where 7 = ( 71 , 72 ) = 
( 1 / 2 , 1 / 2 ). Simple computations show that for x(t) = t, with t £ [0,10], we have 


c 




(10 -t) 1 -^) 
2 T (2 - a(t)) + 2 T (2 — (3(t)) ' 


For x(t) = t the functional reduces to 


J(x,T) 



( 2 a(t) — 1 ) dt. 
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In order to determine the optimal time T, we have to solve the equation 2a(T) = 1 . 
For example, let a(t) = t 2 / 2. In this case, since J{x,T) > —2/3 for all pairs ( x,T ) 
and ff(x, 1) = —2/3, we conclude that the (global) minimum value of the functional 
is —2/3, obtained for x and T = 1. It is obvious that the two Euler-Lagrange 
equations ([2]) and © are satisfied when x = x, since 

<93L[x]" ,/3 (f) = 0 for all t € [0,10]. 

Using this relation, together with 


m*' p { i) = o, 

the transversality conditions (J4|) are also verified. 
For our last example, consider the functional 


J(x,T) 



2a(t) - 1 + 



2T(2-a(t)) 


(10 -t) 1 -?® 

2T{2-m) 


31 


dt, 


where the remaining assumptions and conditions are as in the previous example. 
For this case, x(t) = t and T = 1 still satisfy the necessary optimality conditions. 
However, it is now not obvious that (x, 1) is a local minimizer to the problem. 
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